The sensitivity of Graphene 'Snap-through' to substrate geometry 
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Abstract 

We study theoretically the deposition of Few Layer Graphene sheets onto a grooved substrate 
incorporating adhesion between substrate and sheet. We develop a model to understand the 
equilibrium of the sheet allowing for partial conformation of sheet to substrate. This model gives 
new insight into recent observations of 'snap-through' from flat to conforming states and emphasizes 
the crucial role of substrate shape in determining the nature of this transition. Our analytical 
results are consistent with numerical simulations using a van der Waals-like interaction . Finally 
we propose a novel substrate shape that should exhibit a continuous, rather than 'snap-through', 
transition. 
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Recently, considerable research effort has focussed on characterizing the mechanical prop- 
erties of Few Layer Graphene (FLG) sheets - from thefr elasticity and strengthflQ to their 
ultrastrong adhesion to substrates [3|, |4[ . Understanding these properties is important for 
a number of potential applications since they influence the form taken by sheets, which in 



turn can influence its electrical characteristics 



5]. While standard adhesive tests have been 



applied with some success to measure adhesion, it has also been proposed that deposition 
onto a corrugated substrate may provide a simpler assay |(|. For a given substrate geom- 
etry and material properties, we expect that a sheet deposited onto such a substrate will 
adopt one of three configurations (illustrated schematically in figured]). For relatively weak 
adhesion, we expect the sheet to sit above the substrate with very little deflection (the 

nn 

nonconformal[7|, [8| scenario in figured]). For very strong adhesion, we expect the sheet to be 
significantly deflected and to adopt essentially the form of the substrate (the conformal % |8| 
scenario shown in figured]). Previously, it has been assumed that the transition from non- 
conformal to conformal morphologies is sudden, leading to this transition being referred to 
as 'snap-through'. Indeed, such a snap-through transition has recently been observed in 
FLG sheets 9|. However, in principle a third, intermediate, morphology exists, which we 
term 'partially conformal' (see figured]), with the sheet conforming to the substrate over a 
finite portion of its length but not everywhere. In this Letter, we study the transition from 
nonconformal to conformal morphologies theoretically. The question of principal interest 
is whether this transition is sudden (i.e. occurs at a critical adhesive strength) or, rather, 
whether there is a range of adhesive strengths for which a partially conformal morphology 
may be observed. 

We model the FLG sheet as an elastic beam with bending stiffness B and thickness 
h whose position is given by z = w(x). To allow for analytical progress we consider a 
substrate with a single two-dimensional groove of width, /, and depth, 5; the detailed shape 
of the groove is given by z = w s (x). To determine whether the sheet is conformal, we must 
determine the shape of the sheet w(x) and compare this with w s (x). In regions where the 
sheet is conformal with the substrate we have, by definition, that w(x) = w s (x). In regions 
where the sheet is not conformal with the substrate its shape satisfies the beam equation 



where, for simplicity, we neglect the possibility of a tension within the membrane. (The 




FIG. 1: (Color online) Schematic representation of the three possibilities for FLG morphology on 
a grooved substrate. Left: no adhesion and the FLG sheet remains approximately flat, Centre: 
partial adhesion in which the FLG sheet is out of contact with the substrate in the region A^ < 
X < (1 — X)l, Right: complete adhesion with FLG sheet conforming to substrate morphology 

neglect of the tension within the FLG sheet amounts to neglecting a frictional interaction 
with the substrate; the inclusion of such a friction would complicate the analysis and require 
the ad hoc assumption that the sheet be on the point of sliding everywhere.) We assume 
that the groove is symmetric about x = 1/2 and denote the position of the contact points 
between sheet and substrate by x = XI and x = (1 — A)/. The shape of the FLG sheet is 
thus 

{a + a 2 (l/2- x) 2 , \l/2-x\<\l 
(2) 
w s (x), \l/2 - x\> XI 

where the constants do, «2 and the value of A are determined by boundary conditions that 
we will discuss shortly. 

We envisage that the deflection of the free portion of the sheet is caused by the adhesive 
interaction energy per unit area, 7, between the substrate and sheet. In particular, we expect 
that the value of A (which determines the contact points) will be determined to minimize the 
energy of the system. This energy, U, comprises of the bending energy of the sheet (caused 
by its curvature, which drives the sheet to remain nonconformal) and the energy released by 
the sheet coming into contact with the substrate over a portion of its length, which drives 
the sheet towards being conformal. We have that 

nXl nl/2 

U= [w"(x)] 2 dx+ / K'(x)] 2 dx-2A/ 7 , (3) 
Jo J XI 
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where we take the nonconformal state as the ground state of energy. Using the calculus of 



variations, it can be shown 
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11] that the solution w(x) that minimizes the energy ([3]) is 



given by the solution of the beam equation <Q subject to the boundary condition that[20j 



w"(Xl) - <(AZ) = y f 2^/B = V2/£ ec , (4) 



where £ ec = (-B/7) 1 / 2 is the elasto- capillary length [12j . The continuity of sheet displacement 
and slope at the contact point give the constants a and a 2 in as 

a = w s (\l) - (5l 2 {l/2 - A) 2 , 
a 2 = -w' s (\l)/l(l-2X), 

which may be substituted into (j3]) to give a single equation for A for given values of the 
substrate geometry. 

In what follows, it will be useful to rescale vertical dimensions by the depth of the 
substrate groove, S and horizontal ones by its width, I, i.e. W = w/8, X = x/l, etc.. This 
non-dimensionalization introduces (via (Jl|)) the dimensionless strength of adhesion 

T = (l 4 /S 2 ) 7 /B = (R/£ ec ) 2 , (5) 

where R = I 2 / 5 is the typical radius of curvature of the substrate. In physical terms the 
parameter T tells us whether the adhesive energy is strong enough to overcome the bending 
energy penalty resisting the sheet conforming to the substrate. It still remains to be seen, 
however, whether, for a given substrate shape W S (X) the transition between nonconformal (A 
= 0, small T) and conformal (A = 1/2, large T) is smooth or, rather, a discontinuous 'snap- 
through' transition. We couch our study in terms of varying the dimensionless adhesion 
strength T, which may be varied by fixing the bending stiffness B and varying 7 or by 



holding 7 constant and varying B 13[. In the case of FLG sheets the latter approach has 
been achieved experimentally by varying the number of molecular layersjsj]. 

We note that the dimensionless version of the boundary condition eq. (jlj), may be written 
explicitly in terms of the substrate geometry as 

( 2 r)V2 = -W'J{\) - 2W a '(A)/(l - 2A) = W(A). (6) 

For a given substrate shape W S (X) and dimensionless adhesion strength T we therefore have 
a single equation for A (the position of the contact points). Here we shall consider three 



substrate morphologies to illustrate some of the different behaviors that can be observed: 



W„ 



1(1 



w. 



(2) 



+ cos 2ixX) 
sin 4 7rX, 



(7) 



W, 



(3) 



l-(l-ll-^l) 3 , 



for < X < 1 (outside this interval W s = 1 in each case). For each of these substrate shapes 
it is a simple matter to plot the behavior of the RHS of (j6j as a function of A (fig. [2]). We 
see that for a given value of T there are typically 0, 1 or 2 values of A that satisfy ([6J). The 
different behaviors shown in figure [2] influence the nature of the transition from nonconformal 
to conformal morphologies, as we shall see shortly. 




FIG. 2: Main: The function W(A) given by © for the three different substrate shapes considered 
here: (red), (green) and (blue). The intersection of these curves with the horizontal 

black line, (2F) 1 / 2 , illustrates the values of A satisfying Q; these points correspond to maxima 
(dashed curves) and minima (solid curves) of the energy. Inset: the three substrate shapes (colors 
as in main figure). 



The roots of equation (J6j correspond to the extrema of the energy. However, as well 
as a minimum of energy (corresponding to the contact point we would expect to observe 
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experimentally), there may also be other extrema i.e. maxima and inflection points. Eqn ([6]) 
does not contain information as to which of its solutions correspond to minima and which 
are rather maxima or inflection points. In principle, it is possible to obtain this information 
from the variational approach, using the second variation[l^. However, it is also possible 
for global energy minima to occur with A = 0orA=l/2; since these are not extrema they 
do not appear as solutions of flSj) and must be detected by considering the dimensionless 
energy U = U/(B5 2 /I 3 ), with U as in ([3]). We therefore detect the nature of all extrema by 
considering U. 

We shall shortly discuss the different equilibrium states of an FLG sheet for different sub- 
strates using the formulation above and considering the value of A for which U is minimized. 



BE 



151 ]. is to 



However, an alternative approach, which has been adopted in related studies 
model the molecular forces of adhesion directly by means of a medium range attractive and 
short range repulsive van der Waals force. For a thin sheet resting on a nonpolar substrate 



the interaction energy between a sheet molecule and a substrate molecule is given by 16] 



C/ m _ m = -C (r- 6 - Z>- 12 ) , (8) 

with r the distance between the two molecules and C and D material dependent parame- 
ters. To obtain the full interaction energy between the FLG sheet and the substrate, (jHJ) is 
integrated over the semi-infinite substrate and the full thickness of the beam. To simplify 
the resulting expression, we assume that the typical slope of the substrate, 5/1 <C 1. For a 
sheet of thickness H = h/5 and with the distance between the mid-plane of the sheet and 
the surface of the substrate denoted by Y(X) = W(X) — W S (X) the beam equation §I§ is 
modified to become in dimensionless form[? ] 

= {(Y - H/2y 3 - (Y + 

(3[(Y-H/2r 9 -(Y + H/2y 9 ]}. (9) 

The dimensionless constant /3 is related to the equilibrium distance, Yq, between an unde- 
formed sheet and a flat substrate. We define the distance between the bottom surface of 



the sheet and the surface of the substrate in this equi 



ibrium as Y* = Y — H/2. In exper- 



iments presented previously joj], h > 6 nm, y* < 3.3 A JjJ and 5 = 120 nm so that Yq 1, 
V* < ^ « 1. In the limit H > Y* we find that ~ Y* 6 . Furthermore the constant a 
and the adhesive energy T are related by a ~ — (8/3)H 2 T in this limit. (Note that the limit 
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Yq 3> 1 leads to a different simplification of the model, which is appropriate for substrates 
with small scale roughness)?, 8] rather than the larger grooves considered here.) The ordi- 
nary differential equation can be solved subject to free-end boundary conditions using, 
for example, the MATLAB boundary value problem solver bvp4c. The results of these nu- 
merical simulations can be compared with those of the analytical model by examining the 
total energy U. We shall see that the results of this analysis and the semi-analytic approach 
outlined above are generally in good agreement in the limit 1 ^> H ^> Y*. 

Having outlined our analytical and numerical approaches, we now consider the snap- 
through characteristics of the three different substrates given in (j7|). We shall see that 
each of these substrates illustrates a different type of transition; indeed, the transition from 
nonconformal to conformal may be smooth and not a 'snap-through' at all. 

Fully Discontinuous Snap- Through The substrate morphology in (J7J) is of practical 
interest since it closely represents that used experimentally 9J]. Examining the corresponding 
curve for W(A) in fig. [2] we see that only for 27r 4 194.8 < T < 220.2 do partially adhered 
states represent energy minima. From the energy U for this substrate, we find that a global 
minimum exists with A = 1/2 (i.e. conformal morphology) for r > ir 4 ps 97.4 (which is just 
the bending energy of the conformal state). Since this threshold is significantly below that 
at which the small window of partially conformal states exists, we expect the transition to 
the conformal state to be discontinuous, i.e. a 'snap through' occurs. This is confirmed by 
the numerical results using a van der Waals attraction (see fig. [3]) and confirms previous 
assumptions made in the analysis of experimental results 9}. 

(2) 

Partially Conformal States Before 'Snap- Through' The substrate morphology Ws is 
qualitatively similar to Ws l \ albeit with flatter peaks. However, this modification has a 
significant influence on the behavior of W(A) (see fig. [2]). We see that in this case a partially 
conformal state exists as the local energy minimum, provided that V < 529.3. A calculation 
of U shows that this local minimum is the global minimum for V < 257.1 while for V > 257.1 
the conformal state is the global energy minimum. We thus expect that for T < 257.1 the 
sheet will adopt a partially conformal morphology but that for 257.1 < T < 529.3 the 
sheet may adopt either the partially conformal or the fully conformal state or the fully 
adhered state. Which of these states is realized in practice depends on the details of the 
experimental setup, including dynamic considerations. This uncertainty is illustrated by a 
hysteresis loop in fig. [3b, which is also observed in the numerical simulations because of 
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the use of a continuation scheme. Finally, we note that simply determining the adhesion 
energy required to give U = 0, as calculated previously {9], underestimates the value of T at 
snap-through by around 20%. This discrepancy arises because for this substrate geometry 
partially conformal states may have lower energy than the fully conformal state, depending 
on the value of T. This picture is confirmed by numerical simulations (see inset of fig. (3t). 

(3) 

A Smooth Transition The substrate morphology Ws is qualitatively different to those 
of WP and WP since it has a cusp at the midpoint, X = 1/2. The behavior of the quantity 
W(A) is also qualitatively different to that observed for other substrates. In particular, we 
see from fig. [2] that as A increases towards A = 1/2, W(A) diverges. This divergence means 
that for all values of T there is a unique solution of (E]); in other words, the transition from 
nonconformal to conformal morphologies progresses smoothly through partially conformal 
states; no 'snap-through' occurs. Physically, this happens because the curvature of the 
substrate diverges at the cusp and so no finite adhesion energy can overcome the bending 
energy penalty required to adopt a fully conformal morphology. 

In summary, we have presented a new theoretical formulation for the adhesion of FLG 
sheets onto grooved substrates. This formulation provides a simple explanation for why in 
the experiments presented to date a discontinuous 'snap-through' occurs: the curvature of 
the peaks of the substrate is so great that the bending energy penalty that has to be paid to 
conform partially to the substrate is too large. As a result the sheet can only conform fully 
and even then only once the adhesion energy is large enough. Our new formulation highlights 
the crucial role played by the groove geometry in determining the nature of the transition; 
we have presented a substrate shape for which the sheet conforms partially to the substrate 
before snapping-through discontinuously and another where the discontinuous nature of the 
'snap-through' transition disappears entirely to be replaced by a smooth transition from 
nonconformal to conformal morphologies. Furthermore, we have illustrated that in some 
cases the adhesion strength at which snap-through occurs may not be determined by setting 
U = since partially conformed states may be energetically favourable. These findings are 
supported by numerical simulations of a more traditional van der Waals-type model, which 
converges rapidly to the analytic result in the limit Y* <C H <C 1. 

(3) 

Finally, we emphasize that we have proposed a substrate shape, Ws , for which the 
discontinuous 'snap-through' transition is replaced by a smooth family of partially conformal 
morphologies. This development may allow for the determination of the strength of adhesion 



S 



from a single experiment, since the position of the contact points encodes information about 
the strength of adhesion — it is not necessary to do a whole series of experiments with 
different thickness FLG sheets as performed previously While the shape wP may seem 
difficult to fabricate, we note that a qualitatively similar shape is frequently seen at grain 
boundaries 
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19| . The continuous adhesion transition of this groove shape may also be 



of interest for other experiments on thin layer materials; for example, the high curvature 
seen near the cusp may give rise to plastic deformations making this a simple system within 
which to study the plasticity of FLG sheets. 

This publication is based on work supported in part by Award No. KUK-C1-013-04, made 
by King Abdullah University of Science and Technology (KAUST). TJWW is supported by 
EPSRC. 
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FIG. 3: (Color online) Energy U of an FLG sheet above various grooved substrates as a function 
of adhesion strength V. Results show the predictions based on the analytic model for the global 
energy minimum (solid black curves) and local energy minimum (dashed curves) along with the 
numerical results for the van der Waals-like interaction governed by ([9]) with Y* = 10~ 3 (blue 
curves) and Y* = 10~ 4 (red curves), both with H = 0.05. (a) For a 'snap-through' transition 
is observed analytically (see discontinuity in slope of black curve at V ~ 97.4) and recovered 
in numerical simulations as Y* — > 0. (b) For W s a smooth transition from nonconformal to 

(2) 

conformal is observed with both approaches, (c) For Ws partially conformal states are observed 
before 'snap-through' to the conformal state. Note the onset of a partially conformal morphology 
(A), earliest snap-through (B) and final snap- ftirough (C) giving rise to hysteresis loop (arrows). 
Inset: numerical results for U in this case. 



